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Abstract 



Aims. Recently, the study of globular cluster (GC) color-magnitude diagrams (CMDs) has shown that some of them harbor multiple 
populations with different chemical compositions and/or ages. In the first case, the most common candidate is a spread in the initial 
helium abundance, but this quantity is difficult to determine spectroscopically due to the fact that helium absorption lines are not 
present in cooler stars, whereas for hotter GC stars gravitational settling of helium becomes important. As a consequence, indirect 
methods to determine the initial helium abundance among populations are necessary. For that reason, in this series of papers, we 
investigate the effects of a helium enrichment in populations covering the range of GC metallicities. 

Methods. In this first paper, we present the theoretical evolutionary tracks, isochrones, and zero-age horizontal branch (ZAHB) 
loci calculated with the Princeton-Goddard-PUC (PGPUC) stellar evolutionary code, which has been updated with the most recent 
input physics and compared with other theoretical databases. The chemical composition grid covers 9 metallicities ranging from 
2=1.60x10"* to 1.57x10"^ (-2.25<[Fe/H]<-0.25), 7 helium abundances from Y=0.230 to 0.370, and an alpha-element enhancement 
of [Q'/Fe]=0.3. 

Results. The effects of different helium abundances that can be observed in isochrones are: splits in the main sequence (MS), differ- 
ences in the luminosity (L) and effective temperature (Fefj) of the turn off point, splits in the sub giant branch being more prominent 
for lower ages or higher metallicities, splits in the lower red giant branch (RGB) being more prominent for higher ages or higher 
metallicities, differences in L of the RGB bump (with small changes in Teu), and differences in L at the RGB tip. At the ZAHB, when 
Y is increased there is an increase (decrease) of L for low (high) T^s, which is affected in different degrees depending on the age of 
the GC being studied. Finally, the ZAHB morphology distribution depending on the age explains how for higher GC metallicities a 
population with higher helium abundance could be hidden at the red ZAHB locus. 

Key words, stars: evolution; globular clusters: general 



1. Introduction 

Globular clusters (GCs) are objects orbiting around the Milky 
Way and other galaxies and containing between ~10'* and ~10^ 
stars. Nowadays, ~160 GCs are known in the MiUcy Way, with 
metallicities [M/H ] ^ Ior(Z/X) - logiZ/X)^ between —2.3 
and —0.25 (Harrisll 19961 revision 2010) and with ages r anging 
from ~ 6.5 to -13 Gvr dAparicio et alj|20 01; Chabove' 



Marm-Franch et all 120091; iDotter et aLfl20ia: iMoni Bidin et al 



ical findings show clearly that it cannot be the only parame- 
ter at work. In fact, for many years we have been well aware 
of several GCs characterized by quite different HB morpholo- 
gies in spite of their quite similar metallicity. For this reason, 
several second parameters have been invoked to solve this co- 
nundrum, where, of course, the most d irect, and also the first 
invoked, second paramete r was the age dSandage & WallersteinI 
119601: lRood&Ibenlll968h . This problem has been widely de- 



201 lal) . 

For several decades these objects have been considered as 
an excellent laboratory for the study of simple stellar pop- 
ulations (SSPs) because it was believed that all stars inside 
them were formed from a homogeneous cloud at the same 
time. However, it was recognized early that their horizontal 
branch (HB) morphologies could not be simply parametrized. 
While the metallicity o f a GC is the first paramete r defin- 
ing its HB morphology dSandage & WallersteinI [T960l) . empir- 



sus (e.g. ,Catelan & de Freitas Pacheco 


1993; Leeetal. 


1994 


Stetson et alJll996; Saraiedini et al. 


1997; Ferraro et al. 


1991 


CatelanI 2009alt iDotter et al.,,20 1 Q; ,Gratton et al.,,20 1 a and ref- 



erences therein). 

Ano ther early invoked candidate was the initial helium abun- 
dan ce (ivan den Berghl I1965L 119671; ISandage & WildevI 119671; 
Har twicklll9'68ir but the major difficulty with this element is 



Send offprint requests to: A. A. R. Valcarce 



its extremely difficult measurement, which requires high tem- 
peratures in order to be detected in stars. Moreover, blue HB 
stars hotter than ~ 1 1 500 K as well as stars on the so-called 
extreme HB suffer from metal levitation and helium settling 
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(e.g., iGrundahl etaPUQgOl: Behr 2003|), leaving only a narrow 
band between ~11 500 and ~8 000 K for detecting the initial he- 
lium abunda nce in HB stars. Even so, its measurement remains 
difficult (see IVillanova et aT]l2009l l2012h . Even though helium 
was one of the first candidates, this explanation of the second- 
parameter effect be came less popular aft er the Milky Way for- 
mation scenario of ISearle & ZinnI (Il978h . where the age is as- 
sumed to be the second parameter. In addition, differences in 
helium among GCs also seemed less likely, due to the increas- 
ing evidence that GCs have similar helium abundances (e.g., 
ISalaris et alj |2004). 

However, after the obser vational evidence of multiple main 
sequences (MSs) in GCs dNorrisI 12004 [Piotto et al. 2005; 
iD'Antona et alJl2005l:lPiotto et alj|2007h . the hypothesis that he- 
lium is the second parameter has regained strength, but in this 
case the difference in the initial helium abundance is mostly 
among stars inside the same GC (ID'Antona et al.ll2002i 120051: 
iNorrisI 12004 iLee et alJ 120051: iPiotto et al.ll2005L l2007l among 
several others). In fact, it seems that the combination of age as 
the second parameter and a helium spread among stars as the 
"third parameter" is a promising attempt (Gratton et al. 2010). 

It is worth noting that, although the existence of a cor- 
relation between the chemical abundance anomalies (such as 
light-element anticorrelation s) and the HB morphology has been 
suggested for a long time (Norris 1981: Smith & Norris"1983^, 
[Catelan & de Frei tas Pacheco 1995), only recently has direct ev- 
iden ce of this corr elation be en obtained, as in the case of the GC 
M4 jMarino et aLlbol Ic i Villanova et alj|2012h . 

As can be noted, the hypothesis that GCs are SSPs is being 
ruled out, due to the chemical inhomogeneities among stars of 
the same GC that have been (directly or indirectly) detected and 
associated to different star formation episodes in the same GC 
(iDecressin et al.''2007';'D'Ercole et al.'2008':'Carretta et al."2010^, 
[Conroy & Spergel 201 1: Bekki 201 1; Valcaixe & Catelan 201 1, 
and references therein). 

In addition to the presence of multiple populations in GCs, 
another development makes the availability of evolutionary 
tracks for differe nt hehum abundances highly n eces sary: as 
pointed out bv Catelan & de Freitas Pachecd (Il996|) and Catelan 
(2009b), and also emphasize d re cently by lNataf& Gould 
(i2012.) . iNataf & Udalskil (l20Tl . and lNataf et alJ (1201 lallbh . steT 
lar populations with different enrichment histories may have dif- 
ferent helium enrichment "laws" (i.e., different Y-Z relations). 

In this article, we show the effects of the helium enrich- 
ment in isochrones and zero-age HB (ZAHB) loci in the theo- 
retical Hertzsprung-Russell (HR) diagram. First in section|2] we 
present the Princeton-Goddard-PUC (PGPUC) stellar evolution 
code (SEC) which we used to create a set of evolutionary tracks, 
isochrones and ZAHB loci for different metallicities and helium 
abundances (section[3]). Then, we show the effects of the helium 
enrichment (section]?]) and end with the conclusions (section]5][3. 



2. PGPUC stellar evolution code 

The PGPUC stellar evolution code (SEC) is an updated 
version of the code created by M . Schwarzschi ld and R. 
Harm dSchwarzschild & HarnJ 119651: ] Harm & Sch warzschild 
19661) . and th e n high ly modifi e d by A. V. Sweiga rt dSweigarti 



1971 



- -, ; d by , 

1973L 11994 [1997a^ ISwei gart & Dema rqud )1972t 
ISweig; 



Sweigart & Gross! [1974 II 9761: ISweigart & Catelad 1 1 9981 and 



' All of the models described in this paper can be downloaded from 
the "PGPUC Online" dedicated website, as described in Appendix IaI 



references therein). It is focused mainly on the study of low- 
mass stars, defined as stars that develop a degenerate helium 
core before the onset of helium burning. Since almost all the 
stars living in present-day GCs can be considered as low-mass 
stars, the PGPUC SEC is ideal to be used as the basis to create 
isochrones and ZAHBs for studying the GC properties. 

In the following subsections we first explain some of the up- 
dates to the input physics included in the PGPUC SEC and then 
provide a comparison with other codes. Unless otherwise noted, 
all physical quantities are given in cgs units. 

2A. Opacities 

For high-temperature opacities the original version of our SEC 
(hereafter called PG) used t he Rosseland m ean opacity tables 
calculate d by the OPAL team (ll glesias et al.ll987) . as distributed 
in 1992 (iRogers & Igles ias 1992, OPAL92), where twelve met- 
als were considered. In PGPUC, we used the most recent version 
of the OPAL tables (Iglesias & Rogers 1996, OPAL96), where 
seven new metals are considered, and some errors in the original 
OPAL code were fixed. 

In the PG SEC, cond uctive opacity was c alculated using the 
iHubbard & Lampd (IT969 . for logp < 6) and lCanutol (11970, for 
logp > 6) tables. However, these tables cover a limited mix- 
ture of elements, and their approximations are not completely 
valid when the matter is in a state between gas and liquid, 
which happens in the interior of low -mass stars in the RGB 
phase (Catelan et al. 1996: Catelan 2009b). In the PGPUC SEC, 
we used the C assisi et al. (2007, hereafter C07) calculations for 
the conductive opacities, which cover the partial and high elec- 
tron degeneracy regimes. These calculations include electron- 
ion, electron-phonon, and electron-electron scattering. 

When the temperature drops below log T ~ 4.5, opacity 
sources that are not considered in the high-temperature regime 
begin to be important. These sources include: i) negative H 
ion, which is a large source of continuum opacity; ii) atoms 
that are not completely ionized; iii) molecules that are formed 
at low temperatures; and iv) dust grains, which are impor- 
tant at very low temperatures. These opacity sources are im- 
portant for the outer layers of cool stars, and thus should be 
included for a proper description of most of the evolution- 
ary phases of low-mass stars. In PGPU C, we used the opac- 
ity tables calculated by iFerguson et al.l (2005, F05), instead 
of the opacity tables calculated by Bell (1994, unpublished) 
for the PG SEC using the SSG synthetic spectr um program 
(iBell & G ustafsson 1978; Gustafsson & BeH"l979). The F05 ta- 
bles were obtained with a modified version of the stell ar atmo- 
sphere code PHOENIX (lAlexander & Ferguson] 1 19941) . which 
includes an improved equation of state, additional atomic and 
molecular lines, and dust opacity. 

2.2. Boundary Conditions 

When a stellar model is computed, boundary conditions are re- 
quired to solve the problem. In the center of a star the conditions 
are that mass and luminosity are zero, while in the outer layers 
a relation between the temperature and the optical depth (r) is 
required to know the behavior of the temperature, and then to 
obtain the pressure. 

In a simplified treatment, this relation can be obtained ana- 
lytically, as lEddington](ll926h did, where T\t) = It*^{t+1), 
which defines the effective temperature (Tetj) as the temperature 
where a photon emitted from the depth with r = 2/3 has a prob- 
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ability of 6 ^/3 ^ q.S of emerging from the surface. However, 
this relation uses the gray atmosphere approximation, that is not 
a bad approximation when: i) the opacity is dominated by H , 
ii) local thermodynamical equilibrium is achieved (which is not 
appropriate for low-surface-gravity stellar atmospheres, where 
densities are too low), iii) hydrostatic equilibrium is valid (the 
pressure gradient exactly balances gravity, which is not valid for 
hot stars), and iv) there is radiative equilibrium (radiation is the 
main form of energy transport). 

On the other hand, reliable photospheric boundary condi- 
tions can be determined without these approximatio ns from ac- 
tual stellar atmosphere models (e.g., VandenBerg et al. 2008; 
iDotter et"al]|2008,) . but these calculations need a lot of comput- 
ing time, and due to the many parameters involved, tabulated ta- 
bles covering the full evolutionary track of a GC star do not yet 
exist. Instead of using a theoretical relation with many approxi- 
mations, or one that is very time consum ing, semi-analytical for- 
mulae are often used. PG SEC uses the iKrishna Swamvl (Il966h 
T-T relation, whic h was determine d using the limb darkening of 
the Sun. However. ICatelanI (l2007 b has pointed out that this rela- 
tion does not fit correctly the Krishna Swamy data, suggesting 
the following improved fit0: 

T\t) = ^r^4^. (0.55 + 5.46t-7.72t'-^ -1-5. 12t' - I.OSt"), (1) 
which is now used in the PGPUC SEC. 

2.3. Thermonuclear reaction rates 

The thermonuclear reaction rates taken into account in the PG 
SEC are th e reactions related to hydro gen and helium burning, 
as given bv lCaughlan & Fowleii (11988 !) except for the higher rate 
of the '^C(a,7)'^0 reaction advocated by Weaver & Wooslev 
( Il993h . In this update, all the reaction rates were changed 
from t he analytical formulas determined bv iCaughlan & Fowlen 
(Il988h to the latest available ana lytical formulae of the NACRE 
compilation (Ang ulo et al.ll 19991) . where each rate is determined 
with an extrapolation of the measured cross section down to 
astrophysical energies, following theoretical approximations. 
However, for two important reactio ns, namely ' ^C(a, y)"'0 and 
"^N(p,y)'^0, the rates u sed are from lKunz etall ( l2002l K02) and 
iFormicola eTaP (l2004 F04fl respectively. 

2.4. Equation of state 

To determine the structure and evolution of a star, it is required 
to solve the four basic equations of structure. However, an addi- 
tional equation is needed to solve this set of equations, namely 
the so-called equation of state (EOS), which determines the den- 
sity (p) knowing the pressure (P), the temperature (T), and the 
chem ical compositio n. In this work we change d the original 
EOS (ISweigarllll973h to Irwin's EOS (FreeEOS. llrwini r2007fL 
which takes into account non-ideal effects, mainly pressure ion- 
ization. Coulomb effects, and quantum electrodynamic correc- 
tions, which can be important even at non-r elativistic tempera- 
tures and non-degenerate densities (lHeckleil[T994) . The 20 ele- 
ments considered in FreeEOS are H, He, C, N, O, Ne, Na, Mg, 
AI, Si, P, S, CI, A, Ca, Ti, Cr, Mn, Fe, and Ni, while all molecules 

^ Catelan's published formula was misprinted. 
^ The analytical formula was determined bv llmbriani et alJ ( l2005h . 
All FreeEOS documentation is available at 
|http : //f reeeos . sourcef orge ■ net/ 1 



are ignored except for H2 and Hj . Because of the several time- 
consuming non-ideal effects included in FreeEOS, different ap- 
proximations have been implemented to solve the equations, re- 
ducing the CPU time required. The most accurate FreeEOS op- 
tion (EOSl) takes into account 295 ionization states of the 20 
elements, while EOS la considers some less abundant metals as 
fully ionized. EOS2, E0S3, and E0S4 use a faster (but less accu- 
rate) approximation for the pressure ionization effect; however, 
while EOS2 considers the 295 ionization states, EOS3 and EOS4 
consider all metals fully ionized when log T > 6. The difference 
betwee n EOS3 and EOS 4 is that the latter uses the approxima- 
tion of 'P ols et al.l d 1995 1 P95) to calculate the number of parti- 
cles used to determine the Coulomb free energy, and it is exact 
when all elements are neutral or completely ionized. However, 
the P95 algorithm has a discontinuity problem when H2 is dis- 
sociated, which affects mainly stars with masses below 0.1 Mq. 
In PGPUC we use the E0S4 approximation due to its lowest 
consuming time, and the small difference in almost all the pa- 
rameters with the other approximations (see Table[T]l. As can be 
seen, the mass of the helium core (McHe), one of the most impor- 
tant parameters that defines the whole HB phase, only changes 
by 0.002% with respect to the calculation with EOS 1 . 

2.5. Mass Loss 

Historically, t he mass loss rat e has usually been determined by 
the formula o fl eimer i(|l97lR75), which was obtained empir- 
ically. It reads 

M = 4x lO-'V^A^o/y, (2) 

where L», and are the luminosity, radius, and mass of the 
star in solar units, respectively, and rju is a free parameter whose 
value is usually taken as ~0.5 due to uncertainty in the numeri- 
cal coefficient (e.g., in order to fit the observed HB morphology 
in the GCs the numerical coefficient used is 2 x 10"'^). It is 
now well known that this formula does not provide a good de- 
scription of the available mass loss rates for RGB stars, and sev- 
eral alternative prescriptions have been proposed in the literature 
(e.g., Catelan 2000, 2009b, and references therein). In particu- 
lar, iSchroder & Cuntzi (l2005t SC05) have recently suggested a 
new mass loss rate formula that more properly describes the em- 
pirical data and additional physical constraints. They obtained 

M = 8 X lO-'^;?^ ( -I^f i 1 + -^^] MJyr, (3) 
' M, UOOO/T/ \ 4300^./ ^ 

where L», and M, have the same meaning as in eq. T^ff is 
the effective temperature in Kelvin, is the solar surface grav- 
ity, gt is the stellar surface gravity, and 77 is a free parameter that 
we added to allow varying the final stellar mass, even though the 
value 8 x 10 was calibrated to satisfy the HB morphology of 
two GCs (NGC 5904=M5 and NGC 5927), suggesting that the r] 
value added here is typically ~1 (SC05). The two new parame- 
ters (Teff and gt) are related to the mechanical energy flux added 
to the outer layer of the star (depends on Teff) and to the chro- 
mospheric height (depends on g^). This new formula has been 
tested using well-studied nearby stars with measured mass loss 
rates, reproducing them within the eiTorbars (Schroder & Cunt3 
2007). For these reasons, we used the SC05 mass loss rate for- 
mulation, instead of the R75 one. 
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EOS 


CPU 




L 




R 


Age 




Tc 


Pc 




(days) 






(K) 


(lO'^cm) 


(Myr) 


(Mo) 


(lO^K) 


(10^ g/cm') 


1 


12 


1.00342 


2013.1 


3925.7 


6.758 


13162.68 


0.47977 


7.4432 


9.931 


3 


0.336 


1.00384 


2015.3 


3925.2 


6.763 


13162.73 


0.47978 


7.4435 


9.933 


4 


0.156 


1.00348 


2014.8 


3925.4 


6.762 


13162.75 


0.47976 


7.4442 


9.932 



Table 1. Properties of a 0.8 Mq star (Y=0.25, Z=0.001, [a/Fe]=0.0) at the RGB tip, where the complete evolution was computed 
using the labeled FreeEOS approximation (in case of EOS la and EOS2, PGPUC SEC did not converge with these initial properties). 
The other columns are: CPU time used to calculate the whole evolution from the MS to the RGB tip (based on three Intel Xeon 
E5405 2.0 GHz processors, in parallel-processing mode), luminosity via neutrino emission (Ly), total luminosity (L), effective 
temperature (Teff), stellar radius (R), the time spent from the MS to the RGB tip (Age), the mass of the helium core (McHe), central 
temperature ijc), and central density {pc). 




4 6 8 10 

log(p//'-J 

Figure 1. Fractional error 5 = (2fit - 2numKG07)/2numKG()7 for 
temperatures concerning low-mass RGB stars, where Qfit is the 
neutrino emissivity calculated with the H94 (red lines) or KG07 
(blue lines) fitting formulae, and QnumKCO? is the neutrino emis- 
sivity calculated numerically by KG07. Vertical black lines in- 
dicate the logipljie) value in the center of three stars (Y=0.245, 
and the labeled initial mass and metallicity) along their evolu- 
tion through the RGB phase, when the center reaches the labeled 
temperature. In the panel for log T - vertical lines indicate 
the density for the maximum temperature reached in the center 
(log T ~ 7.87). The log(p/jUe) ranges are determined by the lim- 
its of the KG07 numerical values. 



2.6. Neutrino energy loss 

When a photon is propagating in a plasma under sufficiently ex- 
treme conditions (high temperature, high density), the photon 
can be coupled to this plasma. When this happens, the pho- 
ton can decay into a neutrino and an anti-neutrino (emitted 
in opposite directions), which cannot occur when a photon is 
not coupled to the plasma, since energy and momentum cannot 
simultaneously be conse rved under these circumstances (e.g., 
iKantor & Gusakovl 1200 7'. KG07). In low -mass stars, these ex- 
treme conditions are reached in the degenerate core when the 
star is in the RGB and the white dwarf phases. Generally, in 
stellar evolution codes fitting formulae are used to take into ac- 



count the energy lo ss by this e f fect, th e most common ones being 
those p roposed bv lHaftet all (II994 1995, H94) and lltoh et all 
(1 19961 196). However, KG07 have pointed out the ineffective- 
ness of these formulae to reproduce the real energy loss via neu- 
trino emission for extreme values of p/fi^ while, moreover, the 
196 fitting formula is only valid near the maximum of emissiv- 
ity. KG07 obtained this emissivity numerically, suggesting a new 
fitting formula, which solves the problem for extreme values. 

Figure [T] shows the fractional error of the H94 and KG07 
fits, compared to the numerical neutrino emissivity from KG07. 
It can be seen that the H94 formula reproduces the numerical re- 
sults better in the relevant conditions for low-mass stars (vertical 
black lines). For this reason, we continue using the H94 recipe. 
However, if calculations for extreme values are required, we rec- 
ommend the use of an interpolation of their numerical tables. 

2. 7. Comparison witli other codes 

In this subsection we compare the PGPUC SEC with some 
SECs widely used in the study of lo w-mass stars and 
GCs. These SECs are: Victoria-Regina ( VandenBera et 
l2006i hereafter VIC-REG), BaSTI (Pietrinferni et al. 200^ 
GARSTEC (Weiss & Schlattl 2008), STAREVOL ( Siessll2006l) . 
and PADOVA (.Bertelli et al.. ,2008.) . Table |2] shows the input 
physics used in these codes, where it can be seen that BaSTI and 
GARSTEC include nearly the same physics used in PGPUC. 

In Figure |2] we compare some stellar (internal and exter- 
nal) properties of the SECs that have these data availably 
The star has a mass of O.8M0, an initial helium abundance of 
Y = 0.250 (near to th e primordial helium abundance Yp = 
0.245: llzotov et al]l2007 . and references therein), a metal abun- 
dance of Z=0.001, a solar s caled abundance of alpha elements 
([c/Fe] = 0.0) according to iGrevesse & Noelsl (Il993l) . and the 
evolution is without mass loss. These properties are close to 
the mass, helium, and metallicity of turn-off (TO) stars in GCs. 
Moreover, only for this comparison PGPUC uses the NACRE 
reaction rate for ''*N(p,y)'^0, which is the rate used for all codes 
that use NACRE as input. In that figure the good agreement of 
PGPUC with BaSTI and GARSTEC can be observed, the latter 
being the codes using FreeEOS (interpolated tables). Note that 
the remaining different physical inputs could possibly produce 



^ GARSTEC evolutionary tracks were obtained from the Stellar 
Model Challenge project (http://www.mpa-garching.mpg.de/stars/| 
SCC/Models.htmll, while PADOVA and VIC-REG ones were down- 
load from their respective web pages: http://www3.cadc-ccda.hia-iha. 
nrc-cnrc.gc.ca/community/VictoriaReginaModels and http://pleiadi. 
pd.astro.it STAREVOL evolutionary track was obtained via private 
communication with Ana Palacios. BaSTI evolutionary tracks were 
obtained from the official webpage http://albio nes.oa-teramo.inaf.it/l or 
via private communication with Santi Cassisi. 
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G 5 10 15 5 10 15 



Age [Gyrl Age [Gyr] 

PG BaSTI STAREVOL 

PGPUC GARSTEC 

Figure 2. Comparison of some evolutionary properties for a 
0.8 Mq star with Y=0.250 and Z=0.001, as predicted by some 
of the codes presented in Table |2l PG, dotted black lines; 
PGPUC, black lines; BaSTI, red lines; GARSTEC, blue lines; 
STAREVOL, cyan lines. The properties are luminosity (L/Lq), 
effective temperature (Teff), radius (R), central temperature (Tc), 
density (pc), and hydrogen abundance (Xc). Small boxes are 
zoom-in to show the small differences between BaSTI and 
PGPUC results. 



the observed differences; e.g., if the GARS TEC SEC a dopte d 
the conductive opacities from C07 instead of I toh et all (Il984l) . 
its predicted luminosity at the end of the RGB phase would de- 
crease due to the earlier helium flash (C07). The differences with 
the PG and STAREVOL SECs must be produced by the different 
EOS's which lead to a lower central density (pc) and temperature 
(Tc) than the other three codes, and, consequently, a lower hy- 
drogen burning rate. 

The impact of these differences can also be observed in an 
HR diagram, where in the left panel of Fig. [3] we show the evo- 
lutionary tracks calculated with the same codes as in Fig.jS] As 
can be observed, PGPUC has better agreement in the MS with 
BaSTI, GARSTEC, and STAREVOL, but PGPUC produces the 
hottest MS (AFeff ~ 35 K at a solar luminosity). In the sub- 
giant branch (SGB) these differences are reduced, but at the base 
of the RGB STAREVOL begins to diverge from the other three 
codes. At the RGB bump PGPUC is slightly brighter and cooler 
than predicted by BaSTI and GARSTEC (AT^ff ~ 18 K and 
Alog(L/L0) ~ 0.03 with BaSTI). Finally, at the RGB tip the 
PGPUC track has the lower luminosity. Although there are sev- 
eral other differences along the RGB evolution, they tend to be 
very small. 

For comparison with the VIC-REG model (middle panel of 
Fig. [3]), we produced an additional evolutionary track with the 
same properties as used to produce the left panel of the same 
figure, except for the initial helium abundance that in this case 
is Y = 0.237. However, due to the different input physics be- 
tween PGPUC and VIC-REG, the differences along the evolu- 
tion are relatively larg^. In the right panel of Fig. [3] we also 



* During the publica tion of this paper, an upd ated version of VIC- 
REG was published bv lVandenBerg et"al] ( l2012h . which could reduce 
these differences. 



compared with PADOVA evolutionary track, where the proper- 
ties are: M = O.8 M0, Y ^ 0.230 , Z ^ .001, and [a/Fe] = 0.0 
according to Grev esse & Sauva ll ( Il998l) . It reveals some small 
differences similar to those in the left panel, but at the RGB 
bump there is a difference of Alog(L/LQ) -0.15 (PGPUC is 
brighter). Also, at the RGB tip the difference has the same mag- 
nitude, but PGPUC is instead fainter. This could be due to the 
differences in the conductive opacity tables and in the thermonu- 
clear reaction rates, because both changes produced an increase 
of the RGB bump luminosity and a decrease of the RGB tip lu- 
minosity. 

Due to the fact that this paper is focused on He-enhanced 
theoretical models, in Fig.|4]we compared PGPUC and BaSTI 
results for [Fe/H] = -1.62 with Y=0.30 and Y=0.35, 
where PGPUC results are interpolated as is described in 
Appendix [A] (see a lso Po tter et al,. 2008; Be rtelli et al. 20091 
iDi Criscienzo et al.1 bOlOi: iDeirOmodarme et al.l 120121) . In the 
left and middle panels, we show the comparison of evolution- 
ary tracks with masses between 0.5 and I.IMq in intervals of 
AM - O.IMq. As can be seen, there are not large differences 
at the MS except for 0.7 and O.5M0, where PGPUC tracks are 
cooler. At high luminosities it is possible to observe that PGPUC 
results are also cooler and fainter than BaSTI tracks. In the right 
panels are shown the comparisons between the ZAHB loci of 
the two SECs, which are created using a progenitor mass of 
O.IMq. For this comparison we have to move the BaSTI ZAHB 
loci by Alog(L/L0) = -0.025, since the BaSTI's code uses the 
P99 conductive opacity table inducing a brighter ZAHB than the 
C07 ones (see C07). Apart from previous difference in luminos- 
ity, there is still a difference at the red part of the ZAHB, with 
the most important factor likely being the rate of the reaction 
'^N(p,y)'^0 that is adopted in both codes. More specifically, for 
this comparison the PGPUC SEC uses the F0 4 rate, whereas 
BaSTI uses the NACRE one. iPietrinferni et al.l (l20I0h pointed 
out that two ZAHB loci created with these rates have a difference 
in luminosity of the order presented here at log Teif = 3.85, the 
ZAHB determined with the NACRE reaction rate being brighter 
(at this metallicity). As far as the predicted temperatures are con- 
cerned, we do not find any significant differences. We also show 
the ZAHB loci created with a progenitor mass of O.9M0, and this 
case is discussed in more detail in section l4.4.1l Finally, we con- 
clude that PGPUC models are in agreement with BaSTI ones, 
except for ZAHB loci being fainter by Alog(L/L0) = -0.025, 
due to the differences in the adopted input physics. 

3. PGPUC database 

3.1. Evolutionary Tracks 

To study the effects of a helium enrichment in globular clus- 
ters, we created a database of evolutionary tracks covering the 
evolutionary phases from the zero-age main sequence (ZAMS) 
to the tip of the RGB, together with zero-age horizontal branch 
(ZAHB) loci. This database of evolutionary tracks includes stel- 
lar masses between 0.5 and 1 . 1 Mq with AM -0.1 Mq, for 7 he- 
lium abundances, 9 metallicities, an enhancement of the alpha el- 
ements by [a/Fe] = 0.3, and a mass loss rate according to equa- 
tion [3] (with T] = 1.0). The metallicity values were selected to 
obtain the iron abundances [Fe/H] shown in the second column 
of Table |3] when Y - 0.245. Since it is assumed that nearly all 
GCs do not undergo internal chemical enrichment due to super- 
nova explosions, the mass fraction of Fe within a given GC (but 
not [Fe/H]) woul d be constant over the rang e of possible helium 
abundances (e.g- IValcarce & CateIanll20Tll) . For this reason we 
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dSchroder & Cuntz 2005), VBOO ( VandenBer g et al.ll2000l) , W W93 jW eaver & Wooslev 1993h . Chemical compositions: GN93 
(Z/X = 0.0244. ,Grevesse & Noels., 1993,) . GS98 (Z/X = 0.023 l. lGrevesse & Sauval.,1998i) . 






PCPUC (GN98)- 
- PADOVA (GNSIB)- 



Y = 0.230 



3.9 3.8 3.7 3.6 3.9 3.8 3.7 3.6 3.9 3.8 3.7 3.6 

-og (T.„) Log {T.,J Log {?.„) 

Figures. RGB evolutionary track comparison between PGPUC and the labeled codes (see Table |2]i for a 0.8 Mq star with three 
different helium abundances: Y=0.250 ( left), Y=0.237 (middle) a nd Y=0.230 (right). The metallicity is Z=0.001 and the chemical 
composition is solar-scaled according to'Grevesse & Noelsl (Il993h . except when labeled GN98, where the chemical composition is 
solar-scaled according to Grevesse & Sauval ( 1998). Each panel contains an upper left and a lower right panel with zooms to the 
RGB tip and RGB bump regions, respectively, with the same scales for each Y. 



will use the same metallicity value for each set of Y values in our 
analysis. Even though variations of other elements can be asso- 
ciated to He-enrichment and also produce some changes in the 



theoretical evolutionary tracks and ZAHB loci (C-hN-hO and/or 


other heavier elements; 


Milone et al. 2008; Cassisi et al. 


2008; 


Pietrinferni et al.l'2009; 


Ventura et al. 2009; Marin-Franch et al.i 


2009; Alves-Brito et al 


2012; VandenBerg et al. 2012), 


in this 



paper we focus exclusively on He. We leave the calculations with 
non-standard metal ratios for future papers, where we also plan 
to increase the ranges of masses, metallicities, and [(r/Fe] ratios 
covered. 



We adopted the sola r chemical abundances presented by 
iGrevesse & Sauvall dl998i GS98), where Z/X = 0.023 1 + 0.005, 



when calibrating the modelfl Using this solar composition, we 
computed a Standard Solar Model, assuming an age for the Sun 
of 4.6 Gyr. By fitting the present solar luminosity, radius and 
Z/X ratio, we obtained a mixing length parameter a - 1.896, 
an initial solar helium abundance Yq = 0.262 and a global solar 
metalUcity Zq = 0.0167. 



' Although recent anal yses using 3D hy d rodynamical models 
of the solar atmosphere jXsplund et^ 120051 : iLodders etal] 120091 : 
jZ^affau et al. 201 1) have resulted in lower values for Z/X compared to 
the previous GS98 mix, we do not use these lower values because they 
disagree with high-precision helioseismol ogical determinations of the 
solar structure (e.g., Basu & Antia 2004; iSerenelli et al.l 120091 , 120111 : 
iBi et alj201 iklAntia & Basu,201 L and references therein). 
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Figure 4. Comparison between PGPUC and BaSTI results for non-standard chemical compositions. Left panel: Evolutionary tracks 
with Y = 0.30 and Z = 9.2 x 10 Middle panel: Evolutionary tracks with Y = 0.35 and Z = 8.5 x 10 Right panels: ZAHB loci 
for Y = 0.30 and Z = 9.2 X 10 (upper panel) and Y = 0.35 and Z = 8.5 x 10 (lower panel), where for comparison BaSTI ZAHB 
loci ai-e moved A log(L/L0) = -0.025. PGPUC ZAHB loci are shown for progenitor masses of M,- = O.TMq and M,- - 0.9Mq. 
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Table 3. Iron abundances [Fe/H] as a function of the metallicity value Z and the helium abundance Y for an (T-element enhancement 
of [Q'/Fe]=0.3. In boldface type are the selected [Fe/H]'s (when Y = 0.245) to obtain the Z values used to calculate our theoretical 
models. 



3.2. Isochrones 

For studying GCs, if the canonical SSP is accepted, it can 
be assumed that all stars were formed at the same time with 
the same chemical composition. Since higher-mass stars evolve 
more rapidly, one cannot compare the evolutionary tracks for an 
SSP directly with the GC observations. Instead, it is necessary to 
build isochrones which connect the points at the same age along 
the evolutionary tracks with different masses. 

To construct our isochrones, we use the e quivalent evolu - 
tionary points (EEP) technique first proposed bv iPrathe^ (1 1 976h . 
These EEPs represent points with similar properties (which can 
be internal or external) for tracks of different masses. For low- 
mass stars, the EEPs are: 



- EEPl: This point is defined as the beginning of the MS, 
where we take, as a first approximation for low-mass stars, 
the age to be zero (for SECs including the pre-MS evolution, 
this point is defined when the luminosity due to gravitational 
collapse decreases to 0.01 Lq). 

- EEP2: This point is defined when the helium abundance in 
the center reaches a value of 90%, which is near the point 
when low-mass MS stars have reached their maximum ef- 
fective temperature. 

- EEP3: This point is similar to EEP2, but is determined when 
the central hydrogen is completely depleted. 



- EEP4: This is the point where the mass of the convective 
envelope is 25% of the total mass, which characterizes the 
base of the RGB. 

- EEPS: This point is determined when the star's luminosity 
starts to decrease as the H burning shell arrives at the chemi- 
cal composition discontinuity previously produced when the 
envelope convection zone reached its maximum depth. 

- EEP6: This is where the star's luminosity starts to increase 
again, as the stellar structure readjusts after EEP5 to the 
higher H abundance outside the composition discontinuity. 
For the cases where the RGB bump does not exist (M < 
0.7 Mq in Fig.|5]), EEP5 and EEP6 are determined by extrap- 
olating the RGB bump positions for lower-mass stellar tracks 
with well-defined RGB bumps. 

- EEP7: Finally, the last point is defined when log(LHe/i'0) > 
2.0 due to the beginning of the helium flash or when Teft 
starts to increase due to the star's leaving the RGB before 
igniting He in the core. 

The definition for these points was kindly provided by S. 
Yi (2008, private communication), except for EEP3 which was 
added to obtain a better resolution, and EEP7 where the origi- 
nal definition took into account only the luminosity via helium 
burning. Indeed, the original definition of EEP7 was valid only if 
the star could ignite helium, whereas our new de finiti on is valid 
for any star (see also lBergbusch & Va ndenberg |I992h . Figure |5] 
shows an example of some isochrones between 7 and 15 Gyr, 
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Figures. Set of evolutionary tracks (black lines) for stars 
between 1.1 (hottest track) and 0.5 Mq (coolest track) with 
Y=0.245, Z=9.3xl0-4, and [a/FeJ^+O.S, plotted at AM = 
O.IMq intervals. Red points are the EEPs for each track, while 
green lines are isochrones, from left to right, between 7 and 15 
Gyr, in steps of 1 Gyr. Inner panels are zoom-ins to the RGB tip 
(left) and RGB bump (right). 



which were created using PGPUC evolutionary tracks, the tech- 
nique presented before, and the Hermite interpolation algorithm 
for constructing reasonable analytic curves through discrete data 
points presented by Hill ( 1982). 

Finally, we emphasize that the EEP technique can also be ap- 
plied to tracks for different Y, Z, [a/Fe], and any other parameter 
(mixing length parameter, mass loss rates, overshooting, etc.). 
This means that one can always create new evolutionary tracks 
with a parameter inside the original set of tracks, maintaining 
the others constants. For example, instead of interpolating in 
mass (as was used to create an isochrone), one can interpolate 
in Y to create a new evolutionary track for a different helium 
abundance while maintaining fixed the other parameters (ini- 
tial mass, Z, [a/Fe], C-i-N+0 abundance, etc.). This procedure 
can be repeated for each parameter, with the end result that one 
can always create a new evolutionary track (and consequently 
isochrones) with any parameter inside the grid of original param- 
eters. Such a procedure is implemented in the PGPUC Online 
web page (see Appendix lAll. 

3.3. ZAHB loci 

At the moment when the heUum flash and the subsequent sec- 
ondary flashes have completely removed the degeneracy in the 
helium core, the star has arrived at the ZAHB. In this phase, 
the star is burning helium in the core, while in an outer shell it 
keeps burning hydrogen. Because stars in the HB phase spend ~ 
50% of their life without varying their luminosity by a large frac- 
tion (< 5%, based on PGPUC calculations), it is useful to create 
ZAHB loci which represent the lower boundary of the HB loci. 



ISerenelli & WeissI (l2005h have compared different methods 
to construct ZAHB models without computing the complete evo- 
lution for several stars with different masses and mass loss rates. 
They concluded that ZAHB loci could be well-defined if a star 
is evolved without mass loss from the MS to the HB (passing 
through the flash phase), and then stopping the temporal evolu- 
tion of this star at the point when the rate of change in the in- 
ternal thermal energy is smallest. At this point the carbon abun- 
dance in the center is 5%, and the core convection has reached 
the center. This defines the red end of the ZAHB sequence. Then 
mass from the envelope is removed until the minimum envelope 
mass is reached, which corresponds to the removal of almost all 
of the envelope at the hottest point along the ZAHB sequence. 
The ZAHB locus corresponds to the line connecting all these 
ZAHB models, computed for different envelope masses but with 
the same core mass. We wil l there fore adopt the method sug- 
gested by 'Serenelli & WeissI (l2005h when constructing ZAHB 
loci with the PGPUC SEC. 

4. Effects of AY 

As is pointed out in |T] a spread in helium abundance among 
stars has been suggested to explain why some GCs show mul- 
tiple populations in their CMDs, which is difficult to verify be- 
cause the initial abundance of this element can only be measured, 
in such old populations, at high resolution and with high S/N, 
over a small range in temperatures along the HB (IVillanova et alJ 
120091 12012'). In this section, we show the effects of helium en- 
richment upon evolutionary tracks, isochrones and ZAHB loci 
using PGPUC models, which give us some clues of what we 
should observe in GCs color-magnitude diagrams when a he- 
lium enhancement is present. 

4. 1 . Effects of AY on [Fe/H] 

In Table[3]we show how [Fe/H] increases with increasing Y for a 
constant Z. Thus if multiple populations in a GC are only formed 
with enriched helium and without any enrichment in metals (e.g., 
Valcarce & Catelan 201 1) . sta rs with different Y will show dif- 
ferent [Fe/H] values (see also ISweigartj|l997bt iBragadia et al.l 
'2010a'). For example, if in a GC a population exists with a 
primordial initial helium abundance (Ypc =0.245) alongside an- 
other population with Ysg =0.345, and if both populations have 
the same metallicity Z =1.6 x 10"^ (corresponding to [Fe/H] = 
-1.25 when Y - 0.245), then one should observe a spread in the 
iron abundance by A[Fe/H]~0.063. In fact, the analytical for- 
mula for this spread in [Fe/H] is given straightforwardly by 

A[Fe/H] = [Fe/H]FG - [Fe/H]sG = log , (4) 

where the subscripts FG and SG refer to the first and second 
generations, respectively. 

From spectroscopic observations iBragaglia et al.l 
have found a spread of A[Fe/H] =0.0 70 among the RGB stars 
in NGC 2808 (average FFe/HI =- 1 . 15. ICarretta e t al."2009'). im- 
plying a spread in helium of AY~0. 1 from equation |4] for a 
metallicity of Z=2.18 x lO""*. This helium spread is close to the 
value suggested to reproduce the differe nce between the bluest 
and reddest MSs in this particular GC jP' Antona et all 120051: 
iPiotto et alJl2007l) . 

On the other hand, in addition to the already discussed pos- 
sibility of a variation in the CNO abundances, and thus different 
metallicity at fixed [Fe/H] for the enhanced populations, there is 
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Figure 6. Evolutionary tracks for stars with 0.9 Mq from the 
MS to the tip of the RGB. The left panel shows the effects pro- 
duced by different metallicities (from left to right: Z=l .6 x 10""*, 
1.6 X 10-^ 5.03 X lQ-\ and 1.57 x lO^^) with Y=0.245. The ef- 
fects produced by different helium abundances (Y=0.245, 0.295, 
and 0.345) for two metalHcities ([Fe/H]=-2.25, and -0.25) ai-e 
shown in the right panel. In parentheses are labeled the age in 
Gyr at the RGB tip for each evolutionary track. 



a small possibility that second generation stars could have higher 
Y and lower Z than the first generation, if there was capture 
of metal-poor material from outside of the GC which was then 
mixed with material processed inside the GC (IConrov & Spergel 
I2011h . 



4.2. Effects of AY on evolutionary tracks 

Before comparing isochrones and ZAHB loci, we show the ef- 
fects of changing the metallicity and the helium abundance on 
the evolutionary tracks for stars with 0.9 Mq (see Fig.|6l). Even 
though many of these effects have been known si nce the first evo- 
lutionary codes appeared (Simoda & Iben 1968: lAizenman et alJ 
Il969t llbenlll974t1Simoda&Ibenl ll970). we summarize them 
again using PGPUC stellar models. Although the input physical 
ingredients have been updated over the years, the main effects 
are almost the same. 

As is observed in the left panel of Fig. |6l stars with lower 
metallicities have higher luminosities and effective temperatures 
from the ZAMS to the tip of the RGB. This occurrence is mainly 
due to the decrease of the radiative opacity when the metallic- 
ity decreases, thus enhancing the outward transport of energy 
from the core to the surface. As a result, metal-poor stars have 
a denser and hotter core, which increases the H-burning nuclear 
rates. This is also the reason for which the core H-burning life- 
time decreases with decreasing metallicity. 

Moreover, the temperature of the RGB decreases faster for 
high metallicities due to the higher opacity and the consequent 
faster increase of the radius. At the RGB bump, two additional 
effects are observed: i) its luminosity decreases for higher metal- 
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Figure?. Relation between age at the RGB tip and the He 
core mass (McHe) reached for stars with different initial masses, 
metallicities, helium abundances and different mass loss rates 
(solid lines: no mass loss, or ?; = 0; dotted line: = 1.0). 
Interpolated lines connect the real values for different initial 
helium abundances Y, from top to bottom: 0.245 (red), 0.295 
(green), and 0.345 (blue). For stars with 77 = 1 .0 and 77 = 0.0, the 
initial masses go, from left to right, from 1 . 1 to 0.5 Mq and from 
0.9 to 0.7 Mq, respectively, in steps of 0.1 Mq. As can be seen 
from the monotonic behavior of the solid lines, the helium core 
mass increases for higher ages (or lower initial masses). In the 
case of stars with mass loss (dotted lines), there is a maximum 
value for McHe because, for the mass loss efficiency assumed 
(t] - 1.0), stars lose a large fraction of their masses at higher 
ages, which decreases the hydrogen-burning rate and leads to a 
premature departure from the RGB; a hot flasher may then take 
place. (See text for further details.) 



licities due to a deeper penetration of the convective envelope 
into the interior during the first dredge-up and ii) its prominence 
also increases for higher metallicities due to the larger disconti- 
nuity in the H abundance, and hence in the radiative opa city, at 
the deepest point reached by the convective en velope (Thomas! 
1967; Iben 1968; Swei gart 1978c ISweigart & Gross 1978; Led 
ll97ilYun&LedlI979h 

Finally, we note that stars with a higher metallicity reach a 
higher luminosity at t he tip of the RGB (Iben & Faulkner 1968t 
iDemarque & M engel 1973). This result is a consequence of two 
competing effects. First, the energy output of the hydrogen- 
burning shell in a more metal-rich star on the RGB is greater 
at each value of the He core mass due to the greater abun- 
dance of the CNO elements which increases the nuclear reac- 
tion rates of the CNO cycle. As a consequence the He core 
then grows and heats up more rapidly, thereby causing He igni- 
tion to occur at a lo wer He core mass (ISweigart & GrossllI978l; 
ICavallo et al.lll998l) . However, the greater hydrogen-shell lumi- 
nosity more than compensates for the smaller core mass at he- 
lium ignition, leading to the increased RGB tip luminosity with 
metallicity shown in Figure |6] 
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Figure 8. Isochrones for (from upper left to lower right) 7.5, 10.0, 12.5, and 15.0 Gyr with metallicities, from left to right, Z\-l .60 X 
10 Z2= 1 .60 X 10"^^, and Z3= 1 .57 x 10 ~ in the log L-log Tgif plane. Different colors represent different initial helium abundances: 
Y=0.245 (red), 0.270 (yellow), 0.295 (green), 0.320 (cyan), 0.345 (blue), and 0.370 (black). 



In the case of different helium abundances (right panel of 
Fig.|6|, helium -rich ZAMS star s have higher luminosity and 7^ 
(Demaraue 1967; Iben & Faulkner 1968; Simoda & Iben 1968, 
1970; Demarque etaL ,1971; Hartwick & VandenBerg. ,1973; 
Sweigart & Gross! i978h . To maintain adequate pressure, i.e., 
hydrostatic equilibrium, a helium-rich star must have a higher 
temperature in the center to compensate for the fewer particles 
which have to support the same mass, with the end result that 
these stars are more compact, which increases the T^ff. In addi- 
tion, the opacity is also lower, and then the greater energy being 
produced in the center can more easily flow outward when the 
helium abundance is increased. 



Even though stars with different Y arrive on the RGB with 
different luminosities, the slope of the track along the RGB in 
a log L-log Teff diagram is very similar. However, there are an 
increase of the luminosity and a decrease of the extension of 
the RGB bump when Y is increased, due to the fact that the 
convective envelope in a helium-rich star does not penetrate as 
deeply during the first dredge-up and therefore the composition 
discontinuity produced at the deepest penetration point is then 
both smaller and further out in mass. 

Finally, the luminosity of the RGB tip decreases when Y is 
increased because the flash occurs at a smaller core mass (see 
Fig. IT). This is a consequence of the faster growth in the core 
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mass and hence faster increase in the maximum core tempera- 
ture. 

4.3. Effects of AY on isochrones 

Figure [8] shows isochrones with ages between 7.5 and 15.0 Gyr, 
for several Z values and helium abundances from Y=0.245 to 
0.370, whereas Figure |9] shows in more detail the RGB bump. 
As can be observed in these figures, there are several effects 
on isochrones when Y is increased whose detailed study would 
be quite time-consuming. For this reason, here we summarize 
them qualitatively, to give us hints of what places in HR dia- 
grams we should focus our attention to, when we tabulate these 
ef fects quantitat i vely; a preliminary discussion can also be found 
in lCatelan et al.l (l2010h . Again, we emphasize that several of the 
following effects were directly or indirectly found in the past 
(e.g . Demarque 1967; Iben & Faulkner 1968; Simoda & Ibera 
1968, 1970;' Demarque etal. 1971; Hartwick & VandenBer? 
1973;.Sweigart & Gross..l978;.Lee.,1978;.Yun & Le&.1979.) . 

- MS: As was pointed out for the case of evolutionary tracks, 
when Y is increased the effects on MS stars with equal 
masses include an increase in both L and Teff. Then, if the 
difference in L between isochrones of the same age is com- 
pared for the same Teff, this is a comparison between stars 
with different masses, where the more massive stars are the 
ones with lower Y. Finally, since low-mass stars and/or stars 
with low He abundance have longer evolutionary times, the 
size of the difference in luminosity between MS stars with 
different Y but the same Z almost does not depend on age in 
the range of ages studied here. 

- TO: When isochrones for the same Z are compared, the TO 
point is affected by the difference in Y and also by the age. 
While an increase in Y drives a reduction of L and an in- 
crease in Teff , an increase in age decreases both L and Teff . As 
can be seen, at the TO point the difference in log Teff between 
isochrones with different Y changes with age at almost the 
same rate independently of the metallicity, but the difference 
in logiL/Lo) changes faster with age for higher metallicities. 

- SGB: Here, the difference produced by an increase in Y 
is observed as a decrease in L. However, this difference is 
more prominent for low ages or high metallicities (see also 
iDi Criscienzo et ari i2010). or in other words, when the dif- 
ference in the initial masses of the stars in this zone is large, 
the difference in the log L-log T^ff loci is also large. 

- Lower RGB: In the lower RGB there are three important 
effects to consider when Y is increased: i) the minimum lu- 
minosity almost does not depend on Y, but depends on Z and 
age; ii) as in the TO point, T^s at the base of the RGB is 
higher for higher helium abundances, while the difference in 
Jett between isochrones with different Y increases for higher 
ages; iii) finally, the absolute value of the lower RGB slope 
is higher for lower helium abundances. 

- RGB bump: The effect of an increase in Y at the RGB bump 
is shown in Fig.|9] where an increase of L and a reduction of 
the extension of the RGB bump can be seen (see Sect l4.2l ). 
However, the maximum T^ff of the RGB bump decreases for 
higher Y, but this effect is more evident for lower metallici- 
ties. An increase in age also increases the difference in L be- 
tween Y values, while the difference in Teff almost does not 
change with age. These effects may be observed as changes 
of the shape of the RGB bump in luminosity functions, where 
for large extensions of the theoretical RGB bump the shape 
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Figure 9. Isochrones for the labeled metallicities at the RGB 
bump. Colors represent the same Y values as in Fig. [8] 

is bimo dal, while for smal l extensions the shape is single- 
peaked (ICassisi et al.ll2002h . 
- Upper RGB: Above the RGB bump, as was pointed out in 
Sect l4.2l the difference in Y does not induce any great dif- 
ference in the RGB slope, the only observable effect being 
a shift in the position of the RGB tip, whose luminosity de- 
creases for higher helium abundances. However, this reduc- 
tion of the RGB tip luminosity (and the luminosity function 
at this position) should be very difficult to observe in color- 
magnitude diagrams of GCs, due to the small number of stars 
near this point. 

4.4. Effects of AY on ZAHBs 

As for the isochrones, in several previous works (e.g., 
Faulkner & Iben l96E; Ibe n & Faulknerl [1961 llben & RoodI 
1970; Sweigart 1987) the effects of different helium abundances 
have been studied in ZAHB loci and HB models for different Z 
values. Here, we summarize the main effects for ZAHB loci. 

Figure [TO] shows changes in the ZAHB locus when the He 
abundance is increased for several metallicities. When Z is con- 
stant, it is observed that for high Y i) cooler HB stars (the more 
massive ones) have high ZAHB luminosities, and ii) hotter HB 
stars (the less massive ones) have lower ZAHB luminosities. 
Cooler HB stars show this effect for the same reason MS stars are 
brighter: high helium abundances decrease the radiative opacity, 
thus stars are more compressed (increasing the internal tempera- 
ture and density), leading to an increase in the hydrogen burning 
rate, while the helium burning rate at the core decreases due to 
the smaller helium core mass (see Fig. |7]i, but the luminosity's 
decrease due to the helium burning is lower than the luminosity's 
increase due to the hydrogen burning, with the end result that the 
total luminosity increases. In the case of hotter HB stars, in view 
of the small He core, together with the fact that hotter HB stars 
have tiny envelopes and thus small hydrogen burning rates, these 
stars have lower luminosities when Y is higher: their luminosity 
comes chiefly from the helium burning rate in the (smaller) core. 
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Figure 10. ZAHB loci for different metallicities (increasing from upper left to bottom right) and initial helium abundances. For each 
Y value, ZAHB loci were constructed using three RGB progenitors with different initial masses (M,), which can be related to an 
age dependence of the ZAHB, which is particularly evident at low metallicities. 



Moreover, the luminosity difference between ZAHBs with 
different Y depends on the metallicity, where this difference be- 
tween cooler HB stars is smaller for lower Z, but the difference 
between hotter HB stars is higher for lower Z. The higher in- 
crease in luminosity for higher metallicity when Y is increased 
is due to the increase of the efficiency of the CNO cycle (due 
to the higher Z) that makes stronger the effect of the increase 
of Y mentioned in the previous paragraph. In the case of hotter 
HB stars, the difference in luminosity is due to the difference in 
the helium core masses, which are greater for lower metallicities 
(Fig. El. 

Although ICaputo & DegFInnocentil ( Il995h have shown that 
the luminosity of the ZAHB at a fixed Tgff depends on age for 
low metallicities, this effect has not been (recently) studied in 
detail for the complete range of T^ff and for different Y. In the 
next section we show the effects of age upon ZAHB loci, which 
must be taken into account for all metallicities. 

4.4.1 . ZAHB dependence on age 

In this section we stress the ZAHB loci dependence on age, 
and how this effect is more important at low metallicity and/or 
high helium abundan ces. Althou gh this effect has been known 
for a long time (e.g., Ilbenlll974h . so far it has been largely ne- 
glected. However, nowadays, due to the high precision of the 
available photometric and spectr oscopic data for HB stars (see, 
e.g., Catelan et al. 2009; Daless andro et al.ll201 itlMoehler et all 
1201 It ,Moni Bidin et al..,201 Ih) . it is increasingly important to 
determine the ZAHB locus for the appropriate He core mass, 
which (as already discussed in Iben's review) depends on the 
metallicity, initial helium abundance and the RGB progenitor 
mass (and thus the age). This is especially true when blue HB, 
extreme blue HB and blue hook stars are studied, and even more 
so when the bluest filters are used. 

In regard to the time required for stars to reach the RGB tip 
(?Tip), it depends on three main factors: i) the initial helium abun- 
dance, where tr,p decreases with higher Y because the H burning 
rate increases; ii) the heavy element abundance, where fxip in- 
creases with higher Z because density and temperature in the 
stellar interior are lower in metal-rich stars, which decreases the 
H burning rate at the MS since the burning rate by the CNO cycle 



is negUgible; iii) the initial mass, since more mass ive stars have 
higher internal temperatures, thus fxi n decreases (llben & Rood! 
IT970I iHartwick & VandenBerg||T973l) . These effects are shown 
in Figure|7] for different [Fe/H], Y, and initial masses. 

As Figure [10] shows, the progenitor mass used to create 
ZAHB loci changes the ZAHB luminosity due to the different 
helium core masses reached by these stars at the RGB tip. The 
helium core mass at the RGB tip is determined by a competi- 
tion between the hydrogen burning luminosity which controls 
the rate at which the helium core grows and heats up and the 
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Table 4. Difference in the bolometric magnitude njboi (four last 
columns) between ZAHB loci with different progenitor masses 
(AM,), for 3 metallicities and 3 helium abundances (first column) 
at the labeled log Tetf . In the third column are listed the difference 
in the ages (in Gyr) at the ZAHB of the progenitor stars used to 
calculate the difference in m\,oi. 
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Figure 11. Relation between age at the RGB tip and the initial 
stellar mass (M,) depending on the initial helium abundance and 
the metallicity. Colors represent different helium abundances. 
Dots, and lines have the same meaning as in Fig. |7] 



cooling rate via thermal conduction and neutrino emission dur- 
ing the RGB evolution. To produce the helium flash, T ~ 10^ K 
and p ~ 10^ g/cm^ are required as condition, which is reached 
when some fraction of the helium core is completely supported 
by degenerate electron pressure. 

As is shown in Fig. |7] the parameters that determine the fi- 
nal helium core mass are the helium abundance, initial mass, 
and amount of mass loss on the RGB - the latter two quan- 
tities depending on the age (for a given chemical composi- 
tion). The sense of the dependence of M^He on Y and ini- 
tial mass is that an increase in either increases the tempera- 
ture at the hydrogen-burning shell, and consequently the he- 
lium flash takes place earlier and with a lower McHe- The de- 
pendence on mass loss becomes relevant only when the total 
RGB mass loss becomes quite substantial, which may lead to 
a decrease in efficiency of the H-burning shell, and hence to a 
He flash that may take place only after the star has akeady left 
the RGB, while the star is on its way to the white dwarf cooling 



curve, or even on the white dwarf cooling curve proper (e.g., 
Castellani & Cast ellanil 119931: iD'Cruz et al.]ll996t iBrown et al " 



2001; MoehleretiDMr 

l]). This effect explains the downward 
variation of McHe with age towards old ages that is seen in Fig. 7, 
particularly for Y > 0.245 and 77 = 1.0. The metallicity also 
plays an important role in the final McHe value, but an increase in 
metallicity induces two effects: i) a decrease in the temperature 
at the hydrogen-burning shell; and ii) an increase of the CNO cy- 
cle efficiency. Then, which one of these effects is more important 
depends on the initial mass of the star, as observed in the change 
of the slope of McHe with Z in Figure |7] The temperature of the 
hydrogen burning shell is determined solely by the mass of the 
helium core. This is the same reason why an RGB star evolves 
to higher luminosities. As in a white dwarf, the higher the mass 
of the helium core, the smaller its radius, but the temperature of 
the outer layer of the core is higher, which increases the temper- 
ature of the hydrogen shell around the core with the consequent 
increase of the hydrogen burning rate (Renzini et al. 1992). 

Depending on the helium core mass, the luminosity of the 
ZAHB will change. FigurefTOlshows this effect in HR diagrams. 




Figure 12. Dependence of the effective temperature of ZAHB 
stars on their mass for three different metallicities: Z[ = 1.60 x 
lQ-\ Z2 = 1.60 X 10-^ andZj = 1.57 x lO^^. 



where ZAHB loci can be misrepresented if the mass used to con- 
struct them does not give the required helium core mass for the 
stellar population being studied. For example, ZAHB loci for 
a GC with [Fe/H]=-1.48 and three populations with Y=0.245, 
0.295, and 0.345 at an age of 12.5 Gyr, must be created from 
RGB progenitor masses of -0.83, -0.75, and -0.68 Mq, re- 
spectively. However, if the age is 7.5 Gyr, the RGB progenitor 
masses for these ZAHBs must be -0.95, -0.87, and -0.79 Mq. 
Then, when ZAHB loci are used together with isochrones, the 
initial mass used to create the ZAHB must be defined by the ini- 
tial mass of the star at the RGB tip of the respective isochrone, 
which depends on Y and Z (see Fig. fTTT i. This means that the 
whole ZAHB depends on the age of the GC. This property is 
more important for high helium abundance and/or lower metal- 
licities 

In Table|4]are shown the errors in bolometric magnitude mboi 
that can be induced if ZAHB loci were created with a progenitor 
mass of 0.9 Mq instead of a progenitor mass determined accord- 
ing to the age of the GC. When GCs with multiple populations 
are studied, the mass of the progenitor mass of the ZAHB has to 
be defined by Fig.[TT]if it is assumed that all the stellar popula- 
tions were born at almost the same time. Otherwise, the mass of 
the progenitor of the ZAHB for the younger population has to be 
defined by Fig.[TT]but with M,-Age relation of the given helium 
abundance displaced in the time that this population required to 
be formed. 
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Figure 13. Effective temperature of the stars arriving to the ZAHB for ages between 4 and 15 Gyr for 3 initial helium abundances 
Y= 0.245 (coolest continuous line), Y= 0.295 (dotted-dashed line) and Y= 0.345 (hottest continuous line), and 3 metallicities: 
1.60 X 10 (left panel), 1.60 x 10"^ (middle panel), and 1.57 x 10"^ (right panel). Dotted lines show the change in the effective 
temperatures when a dispersion of 0.02Mo in the total mass loss is considered. Dashed line shows a age of 13 Gyr. 



4.4.2. ZAHB morphology in multiple populations 

The HB morphology of GCs has been studied for several 
decades, but without a consensus about the second parameter 
(together with the metallicity) to parameterize the HB morphol- 
ogy (see section [1]). Historically, multimodality on the HB has 
been used as an argument in favor of the importance of other sec- 
ond parameter c andidates, besides clus ter ag e (e.g.,LN orris 198 1 ; 
Norris & SmithI (19831 ICrocker e t al.^ Iw, " Roodet al., ,1993i: 
Fusi Pecci et al. 1996, see also|Catelan 2008 for a recent review 
on the subject of gaps and multimodality on the HB). More re- 
cently, with the discovery of multiple populations often extend- 
ing down to the MS, which have largely confirmed the expecta- 
tions based on those earlier HB morphology arguments, HB mul- 
timodality is once again being used as a diagnosti c of the pres- 
ence of multiple populations in GCs (e.g.. lAlonso-Garcfa et alJ 
l20T2h . Here we briefly show that a GC with multiple popula- 
tions could have a unimodal HB morphology, even though these 
populations could have a difference in their initial helium abun- 
dance, as in the po ssible case of 47 Tuc (e.g., Nataf et al. 201 la; 
iMilone et alJl2012) . 

First, Fig. [12] shows the effective temperature of stars when 
they arrive at the ZAHB with a given mass (Mzahb) for three 
metallicities, where there is a range of masses that are concen- 
trated near to the minimum T^ff of the ZAHB, and this range of 
masses is higher for higher metallicities. This figure also shows 
how the relation between T^ff and Mzahb changes with differ- 
ent progenitor masses (upper panel), and different initial helium 
abundances (lower panel). As can be seen, a decrease in the pro- 
genitor mass induces an increase of the minimum T^g of the 
ZAHB, the effect being greater for lower metallicities. This de- 
crease also induces a small increase of the Mzahb for a given 
Teff . When the initial helium abundance is increased the range of 
masses concentrated near the minimum T^^fi of the ZAHB also 
increases, which also induces a decrease of Mzahb for a given 

For estimating the mean Teff of a stellar population one can 
use the final mass of stars at the RGB tip (Mxip) depending on 
the age for a given metallicity and initial helium abundance. 
Assuming that the difference between Mxip and Mzahb is small 



(not a bad approximation considering the short time spent in 
the pre-ZAHB phase), one can estimate Mzahb knowing the 
age of the GC in study, and then estimate the mean T^ff. Figure 
[T3] shows the change in T^ff depending on the age for two co- 
eval initial helium abundances (Y= 0.245 and Y= 0.345) and 
3 different metallicities, where continuous lines show the Teff at 
Mzahb = Mjip, and dotted Unes represent the effect of a pos- 
sible dispersion in the mas s loss rate, which show the Teff at 
Mzahb = Mjip ± 0.0 2Mo (Rood & Crocker 19891: ICatelan et al] 
I2OOII: IValcarce & Catelan 2008). The additional dotted-dashed 
line represents an initial helium abundance of Y= 0.295. As 
can be seen, a 13 Gyr old GC with higher metallicity (right 
panel) can hide populations with a AF = 0.10 with the initial 
helium abundance within a small range of temperatures of the 
red ZAHB, but these population will have a great difference 
in their luminosity. For the same age, metal-intermediate and 
metal-poor GCs (middle and left panel of Fig. [13] respectively) 
cannot hide multiple populations with a large difference in the 
initial helium abundance due to the consequent differences in 
Teff. However, smaller differences in Y can be indistinguishable 
for lower metallicities. Of course, evolutionary effects also play 
a role in the shape of the HB, but for studying these effects HB 
evolutionary tracks and Monte Carlo simulations are required, 
which is not the focus of this paper. 

Finally, we have to recall that here we are using eq. [3] with 
77 = 1.0 to take into account the mass loss during the RGB 
evolution for all stars, irrespective of its chemical composition. 
Unfortunately, there are at present no strong constraints avail- 
able on a possible dependence of 77 on Y. On the other hand, M 
in He-rich populations could in principle be constrained, in view 
of the strong dependence of HB morphology on the He abun- 
dance. Future work along these lines is strongly encouraged. 



5. Conclusions 

In this work, we studied the effects that must be observed in the 
theoretical log(L/L0) -log(reif) plane for isochrones and ZAHB 
loci with different initial helium abundances and fixed metallic- 
ity, which is the most common hypothesis to explain the multiple 
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populations observed in the color-magnitude diagrams in some 
GCs. Our study is based on stellar evolutionary calculations for 
different evolutionary phases in the life of a low-mass star. 

First, we presented the PGPUC SEC that was updated 
with the most recent physical ingredients available in the 
literature, including opac i ties, both radiative a n d conductive 
dlglesias & Rogeri Il996t [Ferguson et al. 2005; Cassisi et al.' 

the rmon uclear reaction rates (Anguloetal. 1999; 
Kunz et all 120021: Formicola et al] 12004 llmbriani et al. 2005); 
equat ion of state (.Irwia .2007.) : mass loss (.Schroder & Cuntzi 
l2005h ; and boundary conditions ( ICatelanl|2007h . Even though 
the neutrino en ergy loss rate was n ot updated, we recommend 
the use of the iKantor & Gusakovl (12007) interpolated tables 
for extreme conditions (high densities), which however are not 
achieved for the stars in our study. 

Comparison between PGPUC SEC and other theoretical 
databases shows similar results. In fact, PGPUC SEC evolution- 
ary tracks and isochrones are in excellent agreement with those 
produced using the PADOVA, GARSTEC and (specially) BaSTI 
codes. 

We create a large set of evolutionary tracks for a wide 
range of metallicities and helium abundances for a GS98 chem- 
ical composition and [afFe]-+Q.3, which are available in the 
PGPUC Online web page. Then theoretical isochrones and 
ZAHB loci were used to show the effects of different helium 
abundances at different evolutionary phases along the HR dia- 
gram, where in isochrones: 

- At the MS, the separation induced by Y increases for high Z. 

- At the SGB, the separation induced by Y increases for high 
Z, but it decreases with the age. 

- At the base of the RGB, the difference in Tgff induced by Y 
increases for high Z and also increase with the age. 

while at the ZAHB loci: 

- When Y increases, the L of the cooler HB stars increases, but 
the L of the hotter HB stars decreases. 

- When the age increases, the L of the hotter HB stars increase 
due to the increase of the McHe- 

- When Z increases, there is an increase in the separation in 
L due to a difference in Y at the cooler HB stars, but the 
luminosity difference induced by the age decreases. 

We emphasized the importan ce of the age in determining the 
ZAHB locus (as pointed out bv ,Caputo & Degrinnocentilll995L 
and references in Sect l4.4.lT l. and also the location of stars when 
they arrive to the ZAHB loci depending on the age. In the next 
papers of this series we will study these effects in different filters 
for comparing with GCs data, and we will also create a new set 
of evolutionary tracks with different CNO ratios. 
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Appendix A: PGPUC Online 

PGPUC Online (http : //www2 . astro . puc . cl/pgpuc j is our 
web page containing all the PGPUC evolutionary tracks and 
ZAHB loci. In this web page the Hermite interpolation algorithm 
for constructing reasonable analytic curves through discrete data 
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points presented by iHilll (Il982h is used, which produces evolu- 
tionary tracks from the MS to the RGB tip for any mass, he- 
lium, and metallicities in the ranges Q.SMq < M < I.IMq, 

0. 230 < Y < 0.370, and 0.00016 < Z < 0.01570, respectively, 
for an [a/Fe] - 0.3 and a SC05 mass loss rate with // = 1.0. 
More especificaly, to interpolate an evolutionary track of a given 
mass Ml, helium abundance Yj, and metallicity Zj within a 
range of masses ({M,), with i=1...7), helium abundances ({Y j}, 
with j=1...7), and metallicities ({Z^}, with k=1...9) we follow the 
next steps: 

- Using the first EEP of each evolutionary track with mass 
M = Ml, helium abundance Y = Yi, and all the 9 metal- 
licities Zk, the first EEP of the track with M = Mi, Y = Yi, 
and Z = Z v is interpolated. This process is followed for each 
EEP of the evolutionary track with M = Mi and Y = Yi 
with the end result that an interpolated evolutionary track 
with M = Ml, Y = Yi, and Z = Zj, is created. 

- Following the previous process for each of the seven masses 
Mi, we have an interpolated set of evolutionary tracks with 
seven masses M,, a helium abundance Y = Yi, and a metal- 
licity Z-Zx- 

- Doing the previous two steps for each helium abundance Yj, 
one has a set of evolutionary tracks with seven masses M,, 
for each of the seven helium abundances Y^, and with a given 
Z = Z,. 

- Now, using the first EEP of each evolutionary track with 
mass M = Ml, metallicity Z^, and all the 7 helium abun- 
dances Y j, the first EEP of the track with Y - Y j, that obvi- 
ously has a Z = Z is interpolated. This process is followed 
for each EEP of the evolutionary track with M = Mi and 
Z = Zx, with the end result that an interpolated evolutionary 
track with M = Mi , Y = Y^^, and Z - Z^ is created. 

- When the previous step is repeated for each mass M, a set of 
evolutionary tracks with 7 masses M,-, Y = Y v, and Z = Z j. is 
obtained. 

- Finally, using the first EEP of each mass M,- one can inter- 
polate the first EEP of the mass Mj. Repeating this process 
for each EEP, one creates the final interpolated evolutionary 
track with M = M^-, Y = Y.„ and Z = Z., (see Fig.rOTi. 

In the case that we have the same set of evolutionary tracks, 
but for different a-element enhancements [c/Fe];, one can fol- 
low the previous procedure for each [a/Fe]/, ending with the 
interpolation of the EEPs for each evolutionary track with M - 
M-c, Y = Yjc, Z = Zv, and [a/Fe]; to create ones with a given 
a-element enhancement. Of course, this procedure can be done 
for each initial property that can be parameterized (e.g., mass 
loss, mixing length parameter, overshooting, C+N+O abun- 
dance, etc.), which is our intention with the PGPUC online web 
page: that the user can select the initial property that he or she 
wants. 

Moreover, since any evolutionary track can be created, in this 
web page it is also possible to create any isochrones within the 
previous ranges of helium and metallicity. However, due to the 
fact that the maximum mass for each chemical composition is 

1 . IM0, the minimum age limit is not constant, with a value equal 
to the required time for the star with 1 AMq of the given chemical 
composition to evolve from the MS to the RGB tip. In the near 
future, we plan to increase the maximum mass value. 

In the case of the ZAHB, they can also be interpolated for 
any chemical composition within the range of chemical compo- 
sitions. This is done using the so-called equivalent ZAHB points 
(EZAHBPs), which are defined using two criteria, depending 
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Figure A.l. Interpolated evolutionary tracks (red lines) from the 
original theoretical set of evolutionary tracks (blue lines). Upper 
panel: metallicity interpolation (metallicity increases from left 
to right) with constant helium abundance Y = 0.245 and mass 
of O.9M0. Middle panel: helium interpolation (helium increases 
from right to left) with constant metallicity Z = 0.0016 and 
masses of 0.9 (hotter) and O.5M0 (cooler). Bottom panel: mass 
interpolation (mass increases from right to left) with constant 
metallicity Z - 0.0016 and helium abundance Y - 0.245. 
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Figure A.2. Intei-polated ZAHB loci (red lines) from the orig- 
inal theoretical set of ZAHB loci (blue lines). Open triangles 
show the locus where the convective envelope dissapears. Upper 
panel: Z interpolation (Z increases from top to bottom) with con- 
stant Y = 0.245 and M,- = 0.9Mq. Middle panel: Y interpolation 
(Y increases from bottom to top at the red ZAHB and inversely 
at the blue ZAHB) with constant Z = 0.0016 and M, = 0.9Mq. 
Bottom panel: M,- (M,- increases from bottom to top at the blue 
ZAHB) with two constant Z (Z = 0.0016 and Z = 0.00016) and 
Y - 0.245. Dots are the initial points for each ZAHB locus, with 
the mass increasing from bottom to top. 



on whether the ZAHB star has a convective envelope or not. 
In Figure IA.2I the open triangles show the points which divide 
stars with and without a convective envelope, cooler stars be- 
ing the ones with convective envelopes. For stars with a convec- 
tive envelope, 34 EZAHBPs are defined depending on the per- 
centage of convective envelope mass, from the coolest point of 
the ZAHB (the star with higher convective envelope mass) until 
the point without a convective envelope. For ZAHB stars with- 
out a convective envelope, 100 EZAHBPs are defined, depend- 
ing on the percentage of envelope mass that ZAHB stars have, 
from the point without a convective envelope (100% of envelope 
mass) until the hottest point at the given ZAHB (0% of envelope 
mass). Using the EZAHBPs it is possible to interpolate directly 
in metallicity and helium abundance (see Fig. lA.2l i. However, the 
method used to interpolate between progenitor masses is more 
complex. First, the EZAHBPs can be used only from the hottest 
ZAHB point until the coolest point of the minimum progenitor 
mass (in this case, 0.7Mq). The cool part of the ZAHB is extrap- 
olated from the ZAHB with the nearest higher progenitor mass 
weighted according to the difference in the stellar parameters be- 
tween the three progenitor masses. As can be seen from Fig. lA. 21 
we note that the interpolated ZAHB loci are quite satisfactory. 
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